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ABSTRACT 


Using the first and the second derivative of flutter velocity 
with respect to the parameters, the velocity hypersurface >s. made 
quadratic. This greatly simplifies the numerical procedure developed 
for determining the values of the design parameters such that a 
specified flutter velocity constraint Is satisfied and the total 
structural mass is near a relative minimum. A search procedure pre- 
sented utilizes two gradient search methods and a gradient projection 
method. The procedure is applied to the design of a box beam, using 
finite-element representation. 

The results of the search procedure applied to a box beam indicate 
that the procedure developed yields substantial design improvement 
satisfying the specified constraint and does converge to near a local 
optimum. 
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CHAPTER I 
INTRODUCTION 

In recent years numerical programming techniques hav.egrown very 
relevant in industrial, scientific and military design. This growth 
has been greatly aided by advances in computer techniques and matrix 
methods of structural analysis. It has become feasible to apply 
optimization techniques to design common but complex structures, 
although most of the development up to this time has been at a 
research level. 

One of the main objectives of opt imum des ign aims at replacing 
intuitive drafting of a structure by analytical methods so that it may 
comply with strength conditions. Thus, designing without the theory 
of optimum design consists of formulation of assumptions and their 
verification by calculations, whereas the optimum design makes it 
possible to determine the exact form directly on the basis of given 
strength conditions. 

The major problem with the above mentioned analysis is that, it 
is time consuming. Most approaches to optimum design are based on 
iterative methods, and the complex structure has to be analyzed at 
each iteration to check the various design requirements before pro- 
ceeding to the next stage of optimum design cycle. This leads us to 
the following two criteria for optimum design. 1) It is essential 
that the method of analysis be a rapid one. 2) The number of 
redesign cycles required to arrive at an optimum should be as small 
as possible. 
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OBJECTIVE 

Keeping the above two criteria In mind, in this report the 
author's aim is to modify an existing optimization scheme to mini- 
mize the mass of a structure subject to flutter velocity constraints, 
as given in (l) . 

For each iterative cycle an assumed quadratic equation for 
velocity hypersurface is used. A relatively simple three bay box 
beam is selected to illustrate the application. The box beam has 
enough design parameters and degrees of freedom to make the problem 
meaningful, and if successfuF can be logically extended to more 
complicated structures. 

In stating the mathematical programming problem it will be 
assumed that the geometry of the aircraft wing has been fixed in 
advance. Let m be the mass of the wing, and Dj , D^, .... 0^ be the 
design parameters, representing cross sectional areas, plate thick- 
nesses, diameters squared, etc., of the wing. It can be seen that 
the mass m is a function of these design parameters, which are selected 
so that the mass m is a linear function of these parameters. Mow n 
side constraints can be placed on these parameters so that they are 
bounded. Letting V denote the required minimum flutter velocity, 
the problem can be expressed as follows: 

Minimize m = m (D., D^, ..., Op) (l“l) 

Subject to 

oj < D. - Oj, i = i, P (1-2) 



if 


and 

V >.V 0-3) 

where v is the flutter velocity corresponding to the design parameters 

D., Dp- The superscripts 1 and u correspond to the lower and 

upper bounds, respectively, v though a function of deisgn parameters, 


is in general nonlinear. 



CHAPTER- I I 


FLUTTER 

Flutter is defined as an aeroelastic self excited vibration in 
which the external source of energy is the air stream. The airstream 
feeds energy into the system by virtue of its position or configuration 
at least as rapidly as it is dissipated by damping (2^ p. 192). 

To clarify the above statement let us place an airplane wing in 
a vacuum. The wing, if disturbed from its equilibrium position, would 
vibrate in its normal modes. These vibration, because of structural 
damping, would slowly die out. If however, the wing is moving through 
air with some constant forward velocity v, and is suddenly disturbed, 
as when a gust of wind strikes the wing, then the subsequent motion may 
be such that the amplitude of vibration will tend to (a) decrease due to 
damping if the air velocity relative to the wing is less than the 
critical speed, (b) remain constant, at the critical air speed (also 
known as the flutter speed of the wing), or (c) increase for a speed 
higher than the critical air speed, which may at times cause':il.istruc- , 
tion of the structure. Because it is easier mathematically to describe 
the aerodynamic loads due to simple harmonic motion, theoretical 
flutter analysis consists of assuming in advance that all dependent 
displacement variables at flutter speed are proportional to where 

(j) is the frequency (real) and i = t and then finding combinations of 
velocity v and frequency w for which this actually occurs. 

A simple example of flutter as given in Bisplinghoff et. a 1 . 

( 3 , p. 528) follows. Consider the case of a rigid, symmetrical airfoil 
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hinged at Its leading edge such that It Is elastically restrained 

from rotating about that edge due to the torsional spring with a 

spring constant equal to ft-lb/rad. This is shown in Figure 1. 

Let I = the moment of inertia of the airfoil about the leading edge 

My = the aerodynamic moment per unit span due to ot(t) 

K = the torsional stiffness of the restraint 
a 

The equation of motion for this slngle-degree-of-freedom system is 


I + K a(t) = M . 

“ dt^ a y 

The flutter condition is solved for by assuming as a solution 


a(t) = e 


imt 


( 2 - 1 ) 


( 2 - 2 ) 


where = constant angular amplitude displacement. 
The natural torsional frequency cu^ is given by 


tsi 


K 

_o 

I 


(2-3) 


Using equations (2-2) and (2-3), and dividing equation (2-1) by 
produces 

r , 

2 


.A 

irpb , 


I + — 


iTpb 

where b = semi -chord, 

p = density of the air, 

hence 


u 

__o 

(0 


M 

■rrpb <i) a 


= 0 


I 


irpb 




- I 


- m =0 

y 


(2-M 
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where m = 


M 

,.Y, . 

T52 

Trpb w a 


represents the dimensionless aerodynamic 


coefficient. For a thin airfoil performing small simple harmonic 


motion in two-dimensional incompressible flow, m^, is given by 


(1, p. 529) 


my = T ^’-a i ‘-h 


where M , L , and L, are complex coefficients which are functions of 
a a h 


reduced frequency k = 


mb 


Equation (2-4) can be separated into real 


and imaginary parts. Thus 


Real {m } = 


TTpb 


and 


Imaginary {m^} = 0. 


(f) 


- 1 


(2-5a) 


(2~5b) 


Thus, the flutter occurs at that value of reduced frequency k which just 
makes the out-of-phase component of the aerodynamic moment vanish, 
provided the corresponding in-phase part is of such magnitude that 
equation (2-5a) yields a real flutter frequency oj (J, p. 529). 

To see how the dynamic instability is caused by the energy added 
by the airstream, consider the work done by the airstream on the simple 
airfoil of Figure 1 as it undergoes simple harmonic motion. Since the 
physical quantities are represented by the real parts of the complex 
representation, 

dW = Real {M };<jReal ^ dt 

y y dt 
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where da = differential increment in a(t) 
and dW = incremental work due to da. 

Therefore,, the total work done during one cycle is 


W = J Real {M }',x.'Real dt. 


Now H =Trpbtoam;m is complex 

y y y 

= Trpb^O)^ (m „ + im ,) a e'^^ 
yR y r o 

where the subscripts R and I denote real and imaginary parts, 


respectively. 


= tan 




,, .^ 2 / I 2 ^ 2 \ iiiit 

M “trpba)( /m„ +m, )e a e 

y \\l y^ y ^ 

Hence 

Real {Hy} = T\(3b ^ ^ Cos (tot + t()) 
= C Cos (cot - Ip) 

where C = "’yl^ 

,, da . io)t 
Now, -rr “ itoa e 
dt o 

= a to (-Sin tot + i Cos tot), 
o 

therefore, 


Real vr = ~ o w Sin tot. 
dt o 
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Substituting equations (2-7) and (2-8) in equation (2-6) 

2tt 

W = -C a w Cos (o)t + i|j) Sin o)t dt 

° ) 

o 


then 


CcL 01 
o 

2 



[Sin (2ait + ^|j) Sin dt. 


W = C a ir Sin ifj 
o 


or 


W = Tia 

o 


Imaginary 



(2-9) 


Since — — is independent of time, the sign of equation (2-9) depends 

lOlt 

e 

upon the aerodynamic coefficients M , L ,and L. . A negative sign 

o( Ob n 

would mean that the airstream is extracting work from the elastic 


system, thus providing aerodynamic damping. A positive sign would 
mean that the airstream is adding energy to the system and would 


cause the airfoil in this simple example to flutter. 



CHAPTER I I I 


WORK BY OTHERS 

Although the technique for analytical prediction of aeroelastic 
phenomena has been available for a long time, to date there have been 
relatively few published works dealing with optimization under aero- 
elastic constraints, since most of the work in optimization has dealt 
with the conventional conditions of strength, stiffness and stability. 

Rudisill and Bhatia (1) have developed a numerical procedure for 
determining the values of the design parameters such that a specified 
flutter velocity constraint is satisfied and the total structural mass 
is near a relative minimum. The search procedure utilized two 
gradient search methods and a gradient projection method. In the 
above procedure substantial design improvement was made but convergence 
to an optimum was not obtained in a reasonable computer execution 
time. Since then Rudisill and Bhatia (/i-) have obtained an analytical 
solution for the second partial derivatives of the eigenvalues of the 
flutter equation along with the equation for finding the second partial 
derivatives of a flutter velocity of an aircraft sturcture with respect 
to the structural parameters. Using these partial derivatives in com- 
puting the step size used in the projected gradient search along a 
constant mass hyperplane has helped in cutting down on the number of 
redesign cycles to arrive at the optimum. 

Cooper (5) in his report, has attempted a direct method of 
solution for the flutter velocity. He also has applied his method 
to the optimization of a cantilevered box beam for minimum mass due 
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to a flutter velocity constraint. His method of solution required 
a plot of the imaginary part of the eigenvalue of the flutter equation 
versus the dependent variable of the flutter equation, i.e. velocity 
divided by the circular frequency. From this plot the crossover 
points on the dependent variable axis are sought. From these points 
the lowest or the critical velocity for which the structure will have 
divergent oscillations may be calculated. He has simplified his pro- 
cedure by fitting *a simple quadratic or cubic equation for the 
required plots, and has obtained good results. 

Siegel (6) within the last year has developed an optimization 
method for accurately and rapidly calculating, through a completely 
automated digital computer program, the minimum weight spanwise dis- 
tribution of an airfoil surface to provide a given required flutter 
speed. This program starts by calculating the flutter speed for a 
configuration having adequate strength. If the calculated flutter 
speed is lower than required an automatically determined increment of 
structural material is added to the spanwise location where the strain 
energy per unit structural volume in the flutter mode is a maximum. 
Using the new structural data the flutter speed is again calculated. 
This process is repeated automatically within the computer until the 
required flutter speed is attained. This method is thus absolutely 
based on the concept that the most efficient structure is one that has 
constant strain energy per unit structural volume In the flutter mode. 
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in this report the method developed by Rudisill and Bhatia (1,^) 
has been extended to include a curve fitting technique for the velocity 
hyperp lanes^ and other minor changes to accommodate it have been madey , 
-t^hus.^ making the process efficient. 



CHAPTER' IV 


DESCRIPTION OF OPTIMIZATION PROCEDURE 

If stiffness distributions are not properly optimized for 
flutter, aircraft structures can become significantly heavier than 
necessary. For the objective formulated in Chapter 1, there are 
numerous search strategies that could be devised to attempt a solution 
to an optimization problem. 

The optimization procedures developed here are basically a 
gradient search method, the greatest rate of improvement of the func- 
tion being found by moving along the gradient. (The derivatives of an 
objective function with respect to each of the n parameters are col- 
lectively called the gradients of the objective function.) 

The present optimization procedure (4) utilizes three well known 
gradient search routines. A simplified graphical illustration of their 
behavior is shown in Figure 3- In the following three sections the 
searches are individually described; later on in this chapter the curve 
fitting method for the velocity hyperplane is given. 

Before we go into the mathematical derivation, a few important 
assumptions are stated. 

(i) the search is to be conducted in the P dimensional space of 

design parameters Dj , D 2 » 0^. The column matrix (D(p)} defines a 

design point p in this space which corresponds to a particular 
structure configuration such that D. > 0, i = 1, P, 

(ii) the flutter velocity is uniquely defined at any point {D}. 
Thus, V which is assumed to be a continuous function of {D} is a 
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scalar point function. The total mass m is also a continuous 
scalar function of {D}. 


Let p and p be two neighboring points defined by {D(p)} and 
{D(p )}, respectively. The scalar distance As between these points is 


As = |{D(p)} - {0(p’')}j = / E (AD )' 

i = l ' 


(A-1) 


It is now assumed that the limits of the difference quotients of v 
and m as As tends to zero, exist. Thus 


Limit 

v(p*) - 

v(p) 

- ^ 

As -»■ 0 

As 


ds ’ 

Limit 

m(p") - 

m(p) 

dm 

As ^ 0 

As 


ds 


Equation (4-2) defines the directional derivatives of v and m. 


VELOCITY GRADIENT SEARCH 


This routine Is employed in order to increase the velocity. In 
this routine the search moves perpendicular to the velocity contours, 
i.e. maximum increase in velocity. The desired velocity is reached In 
one or more steps in an iterative fashion. 

The maximum rate of change of velocity is given by the normal 
derivative and is equal to the absolute value of the gradient vector of 
the velocity. It may be expressed as 


dv 

dn 




(4-3) 
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A new velocity along a gradient curve corresponding to small design 
parameter changes may be approximated by the expression 


V = V + An 
dn 




where V and v are the new and old velocities respectively, 

and An = As along the normal vector to constant velocity hypersurface 

at {D}. The direction cosines of the gradient are 


3v 

3D. 


1 . = 

I 


dv 

dn 


AD. 

I 

An 


. i = 1,P. 


{h-S) 


Therefore, 


3v 

3D. 

AD. 

I 

dn 


An 


, i = 1,P. 


(4-6) 


From equation (4-4) 


(V - v) 

An = 3 . 

dv 

dn 

Substituting equation (4-3) and (4-7) in equation (4-6) 

3v 


(4-7) 


(V - v) 


3D. 


AO. - 


j=1 ^ J 


I 5= 1 P 

9 I — * 9 ' * 


(4-8) 


Thus the new value of design parameters D. , i = 1, P, are given by 


D.“ - D. + 


<V- 1^. 


j-i 


, i = i.p. 


(4-9) 


r 
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Equation (4-9) will yield approximate values of the design parameters 
corresponding to V. 

Since the velocity is a nonlinear function of the design parameters, 
the true velocity corresponding to the set of design parameters D will 
not be equal to V. An iterative procedure may be used to determine the 
set of design parameters corresponding to the desired velocity. The 
search was programmed such that D., i = 1, P, in equation (4-8) is 
either positive or’ zero. 

MASS GRADIENT SEARCH 


This routine is used in order to reduce the velocity. in this 
routine the search moves perpendicular to the constant mass hyperplanes, 
i.e. the direction of the maximum rate of decrease in mass. This 
process is repeated until the velocity is less than or equal to the 
desired velocity. 

The total differential of the velocity may be expressed as 

dv = E <iD,. (4-10) 

The normal derivative of the total mass m may be expressed by the 


equation 


dm ^ / ' X 2 

'"7j=i 


(4-11) 


and the direction cosines of the mass gradient vector are given by 


the relation 


C 
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1 . = 


3m 

3D. 

i_ 

dm 

dn 


AD. 

i_ 

An 


. i - 1,P 


(^- 12 ) 


where An = As along the normal vector to constant mass hypersurface at 
{D}, Approximating the total differential dv by increment v and 
substituting AD. from equation (^-12) in equation (^-lO) yields 


Av 


„ / 3v 3m V . 
.^<3o:3Dr> "" 
J^j ! ! 

dm 

dn 


or 


An 


^ / 3m ^ 2 

1 = 1 I 

' 3v 3m 

. , 3D. 3D.' 

J = 1 J J 

Substituting equation (4-13) into equation (4-12) yields 

3m 


(A-13) 


3D. 


Av 


AD. 


I 


(4-14) 


3v 3m 


. I 3D. 3D. 
J = 1 J J 


The new set of design parameters {D } can again be computed from 

D.'' = D. + AD. 

I 1 1 

It should be noted that Av = V - v. 

The mass gradient search is also an iterative procedure like the 
velocity gradient search. The search was programmed in such a way that 


AD., i - 1, P, in equation would be either negative or zero. 
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GRADIENT PROJECTION SEARCH 


This routine is employed in order to reach a relative maximum 
of the velocity while the mass is held constant. The parameters are 
varied such that the search proceeds tangent to a constant mass hyper- 
surface In the direction of the maximum rate of increase of velocity. 

The gradient projection method allows variations of design para- 
meters to be taken so as to satisfy the behavior constraints at all 
times. The search for an optimum cannot therefore proceed in the 
steepest ascent direction, but must always be restricted to movements, 
satisfying the constraints. Here the velocity v wl 1 1 be maximized 
while the total mass m is held constant. 

With the aid of the second derivatives of the flutter velocity 
with respect to the design parameters, and the projected gradient 
search as given in ["y] , the step size s could be expressed as 


s 


P 

Z 

j = l 


P 

z 

Jrl 


9v 
9D. 
± 


P 

Z 

k=l 


9^v 


9D. 9D 
J 



dPj dD|^ 
ds ds 


( 4 - 15 ) 


dD. 

where are the direction cosines corresponding to the direction of 
the maximum rate of Increase of the flutter velocity along a constant 
mass hypersurface. These direction cosines [ 7 ] are, 



+ X 
^9D. 1 

J 


9m . 

9D.' 

J 


dD, 

ds 


( 4 - 16 ) 
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wh 


ere 2X^ and are given as 


1 - V 3v . " / 9m ^2 

, 3D. 3D. ' r. Sd.'' 
J=1 J J J=^ J 


(4-17) 


«o = 1^1 } 

f- J = 1 J J J ^ 


1/2 


(4-18) 


and m is the total mass of the structure. New parameters may be 
computed from the relation 


Vf 

D. 

J 


“j* ‘1^ 


J 


(4-19) 


If the total mass, is formulated as a linear function of the 
design parameters then the search will always be along a constant 
mass hyperplane except those times when the side constraints for the 
design parameters are encountered. In that case the search might com- 
pletely fail, i.e. velocity decreases instead of increasing. When this 
is encountered a different approach is used which is described later 
on in the computation procedure. 


CURVE FITTING FOR VELOCITY 

In the optimization procedure developed here, we are trying to 
minimize the mass of a structure (aircraft) to satisfy flutter velocity 
requirements. In Its simplest form we will formulate the program for a 
cantilevered box beam, which consists of twelve variable parameters, 
but could be extended to include any finite number of variables. In 
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order to make the procedure systematic and more efficient the velocity 
contours are approximated by a quadratic equation of the form, 


n 


n n 


v=A+ E B.D. + E E C..D.D. 

1 = 1 ' ' 1 = 1 j = l -• 


(i*“20) 


The coefficients A, B., and C.j could now be calculated at the 
known point using the following procedure; Taking the first derivative 
of equation (^-20) with respect to some parameter D|^ 


Sv 


n 


= B. + E C. .0. + E C,. D. 

h hj j . , th I 


n 

E 

i = l 


(^-21) 


and the second derivative with respect to a different parameter D Is 


9^v 


9D. 9D 
h g 


= c. + c . 

hg gh 


/ 


Now assuming above equation yields 


1 


9^v 


hg 2 30. 90 

^ h g 

Substituting this in equation (4-2l) 

n n 

B. = ^ - E C. .0. - E C., D. 
h 9D|^ hj j ih I 


(4-22) 


(4-23) 


(4-23) 


or, 


B. = Ik "2 E C., D. 

h ,h , 


Since C is symmetric, we have 


_ 3v _ ^ 

“^h 90, r. 9C, 90, ^i ' 

h 1 = 1 I h . 


(4-24) 



) 


22 


Again, equation (^»-20) could be written 


A = V 


n n 

I B.D. - Z 


n 


J=1 


E C..D.D. 
1 = 1 j = l 


Substituting values of B's and C's from (4-24) and (4-23) we have 


A = V - Z [ 


3v 


. y 9D 

j=i j 


n 

- E 

. 1=1 


3^v 




3D 




(4-25) 


The first and second derivatives of flutter velocity with respect 
to the design parameters is given in (1 ,^') > the coefficient's 
and A are then calculated with the aid of equations (4-23 through 
4-25). These coefficients are recomputed at the beginning of each 
new design cycle. 


COMPUTER PROCEDURE 

A simplified flow diagram of the optimization procedure is 
shown in Figure 2. and are the specified tolerances used to 
compare the computed velocity v and the computed mass TMASS to the 
desired velocity V and the previous mass TMASSl, respectively, at 
various stages of the optimization procedure. 

Initial parameters are assumed, an initial velocity is computed 

and is fitted to an assumed quadratic surface whose coefficients A, 

B . and C. are evaluated by the method described in the last section, 
h ’ hg 

As shown in Figure 3, if the initial velocity is not within 
tolerance and if v is greater than V, then the gradient mass 
search is executed as shown from a to b. If instead, v is less than 
V, as at point a', then the gradient velocity search a' to h' is 
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executed. The gradient mass search reduces the mass, but may either 
increase or decrease the velocity; the gradient velocity search 
increases the velocity. 

When the velocity is within tolerance e^, then a gradient pro- 
jection search is executed along a constant mass hyperplane for example 
from c to d. If any side constraints are encountered during the 
search, the appropriate design parameters are set equal to their con- 
straint value, and the search deviates from the constant mass hyper- 
plane. if the computed velocity increases from its previous value, 
the search is continued as shown in Figure 3 in an iterative fashion 
until the flutter velocity reaches a maximum for that constant mass 
hyperplane. 

If a large number of side constraints become active, then the 
projected gradient search will deviate greatly from the constant mass 
hyperplane. At that stage the search fails and an alternate procedure 
as shown in Figure A is employed. In this procedure, the gradient 
velocity search and the gradient mass search are alternated in the 
following manner. For the desired velocity V the gradient mass search 
proceeds from a to b Fig. 4, now the desired velocity is changed to 
1.2V and then utilizing the gradient velocity search, a step is made 
from b to c. This procedure is repeated a number of times until the 
change in mass (point O) is within tolerance. The search is terminated 
if the change In the mass for any cycle is less than some small pres- 


cribed number. 



Figure 2 - Simplified Flow Diagram of Optimization Procedure 





















CHAPTER V 


OPTIMIZATION OF A BOX BEAM 

Figure 5 shows a three-bay box beam representing a uniform 
cantilever aircraft wing structure. It Is assumed, as stated in 
Chapter I, that the wing geometry and shape have been fixed in advance. 
Thus the height H, width W and length of each bay L are treated as 
constants, these values could be different for each bay if necessary. 

The design parameters for each bay are defined to be: (I) area of 

longitudinals, (II) front and back web, (ill) top and bottom web 
thickness and (Iv) rib thickness. These parameters were required to 
be of uniform value for each bay. Since there are three bays, the 
total number of variable design parameters Is twelve. Other constant 
parameters needed in the analysis are defined in Figure 5- With the 
design parameters defined as above, the total mass to be minimized is 
a linear function of the design parameters. Solution to the mathemati- 
cal programming problem formulated in Chapter I is found for the box 
beam described above. 

A simplified flow diagram explaining the basic logic of the 
optimization program developed, was explained in Chapter IV (Figure 2). 
Each of the three gradient search routines calls a sub-routine which, 
utilizing the established coefficients of the assumed quadratic, returns 
the computed velocity and the partial derivatives of the velocity with 
respect to the design parameters. In the flow diagram of Figure 2, v 
and V are the computed and the desired velocities, respectively. 

The aerodynamic matrix is formulated from the equation used by 
Smi Ig and Wasserman {S,, p. 398) . 




b => 25 semichord 

a « 2.5", distance of elastic axis from the midchord 


FIGUEJ3 5. Rectangular Three-Bay Box Beam 
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TADDS I IiaTIAL DliSIGK PAR-UiETrlRS 


LPaSi^D: 

J = Kuraber of bajrj 1, 2, 3 
D(J) = Area ox longitudinals for the bay 
D(3+J) = Front and back neb thickness for the bay 
D(6+J) = Top and bottcuu skin thicl<ness for the bay 
D(9+J) = Rib thickness for the bay. 


var;l4BLS 

DESIGN 

?(UIAI'ETSR3 

INITIAL 

VALUES 

CONSTRAINTS 

UPPKR 

m-m 

b(l), Sq. in* 

0.01339 

0.1339 

0.002315 

D(2), 3q. in. 

0.01339 

0.1389 

0.002315 

D(3), Sq. in.. 

0.01339 

0.1389 

0,002315 

D(l.), in. 

0.006667 

0,06667 

0.001111 

D(3), in. 

0. 006667 

0.06667 

0.001111 

D(6), in. 

o. 006667 

0,06667 

o.ooiin 

D(7), in. • 

0.003333 

0.03333 

0 .000555 

D(8), in. 

0.003333 

0.03333 

0.0005h5 

DU')> in. 

0.003333 

0.03333 

0.000555 

D(10), in. 

0.003333 

0.03333 

0.000555 

, D{11), in. 

0.003333 

' 0.03333 

0.000555 

D(12), in. 

0.003333 

0.03333 

0.000555 


origwalpa^ 
OF POOR <50AUlil 








TABLE II. RESULTS OF FLUTTER QPTH‘IIZATIQN 


lEQEND: 

J * Number of bayj 1, 2, 3 

"til 

D(J) = Area of longitudinals for the j bay 

th 

D(3+J) = Front and back web thickness for the j bay 

D^6+J) = Top and bottom skin thickness for 1iie bay 

th 

D(9+J) = Rib thickness for the ^ ' bay 


VARLABI^ 

DESIGN 

PARAIIETERS 

INITIAL 

VALUES 

FINAL 

VALUES 

D(l), 3q. in. 

0,01389 

0.00255U 

0(2), Sq, in. 

0.01389 

0.00231?" 

D(3), Sq. in. 

0.01389 

0.00231?'^ 

D(U), in. 

0,006667 

O.OOIJ4/I 

D(5), in. 

0,00666? 

0.00li4?6 

D(6), in. 

0 .006667 

O.OOllll’'' 

D(7), in. 

0,003333 

0.00055?* 

D(6), in. 

0.003333 

0,0000296 

D(9), in. 

0,003333 

0,00'05988 

D(10), in. 

0.003333 

0.001232 

D(H), in. 

0.003333 

O.C01032 

D(12), in. 

0. >003333 

0.001130 

VELOOm, 



FEET PER SECOND 

670.1 

747.3 

MASS, 



SLUGS 

6.099 

1.098 


-itparameter eaual to the lower constraint 
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The three gradient search procedures developed use partial 
derivatives of total mass and also the first and second partials of 
velocity with respect to the design parameters. When the total mass 
is a linear function of the design parameters, the partial derivatives 
of the total mass are constant and need to be evaluated only once. If 
the total mass is not a linear function of the parameters, they could 
be easily evaluated at each step in the design space. The partial 
derivatives of the flutter velocity are computed as shown in 
reference ,U) . 

The initial variable parameters used for the box beam, and the 
upper and lower side constraints are given in Table I. 

The total mass of the beam and the flutter velocity versus the 
number of redesign cycles are plotted in Figure 6. The desired flutter 
velocity was specified as 800 feet per second. The design was started 
with an initial velocity of 870.1 feet per second and a mass of 
6.093 slug (a. Fig. 6). The velocity was first decreased utilizing 
the gradient mass search, and then increased by the gradient projection 
search, holding mass constant. These two steps were again repeated but 
the gradient projection search failed. An alternate procedure as 
described in Chapter IV (Figure h) was used, and the cycle terminated 
due to no change In mass (b. Fig. 6). The second cycle started with an 
Initial velocity of 388 feet per second. The gradient velocity search 
was used for step one and then the search iterated between the gradient 
projection and the gradient mass search, but did not converge in the 
maximum number of allowed steps. Starting with a velocity of 72^ feet 
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per second and utilizing all three searches it converged to the mass 
of 1.0^7 slugs and a velocity of 770.8 feet per second (d, Fig 6), 
indicating that the program had probably produced nearly a local 
optimum. Three more cycles were carried out with a very small change 
in mass, until finally the search terminated at a velocity of 7^3-7 
feet per second and a mass of 1.098 slugs. Thus a local optimum was 
reached in six cycles, giving the best results so far. 

The flutter velocity was computed by solving the complex eigen- 
value problem expressed in the fundamental flutter equation (l). A 
computer program for solving the complex eigenvalue problem was obtained 
from the National Aeronautics and Space Administration. 

Table 11 shows the results of the optimization study and lists 
the Initial and the final design parameters and the corresponding 
mass and velocity. 

An IBM 370> Model 155 computer at the Clemson University Computer 
Center was used. The total execution time, as recorded by the central 
processing unit (CPU) was 112 seconds for 6 redesign cycles carried out 
by one computer run. Comparing it with a previous search which took 
123 seconds, there is not much saving in time, but this procedure may 
be much more efficient when used with system which has a much larger 
number of parameters and degrees of freedom. 
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CONCLUSIONS 

The quadratic curve fitting for the velocity hypersurface is the 
important feature of this work. The time required for optimizing the 
cantilevered box beam by the method developed here is less than the 
previous method. 

The optimization scheme used is based on three simple gradient 
search techniques, this search have been modified so as to estimate 
step size in case of overshoot. The scheme could be extended to 
include other behavioral constraints such as stresses, displacements, 
and buckling. Since it did result in substantial improvement for a 
system with 12 design variables, the author feels that it is reasonable 
to expect that this scheme would at least work with a similar degree 
of success for a system with larger number of design variables, and in 
fact, might be more efficient. 

Also, there is much room for innovation in the use of search 
schemes which apply other available gradient search methods or their 
combinations to the problems of this nature. 
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